Dipolar gases in quasi one-dimensional geometries 
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We analyze the physics of cold dipolar gases in quasi one-dimensional geometries, showing that 
the confinement-induced scattering resonances produced by the transversal trapping are crucially 
affected by the dipole-dipole interaction. As a consequence, the dipolar interaction may drastically 
change the properties of quasi- fD dipolar condensates, even for situations in which the dipolar 
interaction would be completely overwhelmed by the short-range interactions in a 3D environment. 



Low-dimensional ultra cold gases have recently at- 
tracted a major attention. One and two-dimensional 
gases are created in sufficiently strong optical lattices [H , 
or by means of magnetic wires Low-dimensionality 
leads to a very rich physics, highlighted by the recently 
observed Berezinskii-Kosterlitz-Thouless transition in 2D 
gases [1], the enhanced role of thermal and quantum 
phase fluctuations in elongated gases and the real- 
ization of the Tonks-Girardeau regime of ID bosons . 

The properties of quantum gases are crucially deter- 
mined by the interparticle interactions. Current experi- 
ments typically involve particles at very low energies in- 
teracting via a short-range isotropic potential character- 
ized by an s-wave scattering length a. The latter may 
be modified by means of scattering resonances induced 
by magnetic fields (Feshbach resonance) or by properly 
detuned lasers 0. Interestingly, the scattering proper- 
ties in quasi-lD (also in 2D f^) may be crucially affected 
by the contrained geometry. In particular the transver- 
sal confinement can induce a novel resonance known as 
confinement-induced resonance (CIR) in which the 
transversal trap modes assume similar roles as the open 
and closed channels in usual Feshbach resonances 

Interestingly, new experiments on atoms with large 
magnetic moment cold molecules 11 1 and Rydberg 
atoms [12| have recently opened a fascinating new re- 



search area, namely that of dipolar gases. In these gases, 
the dipole-dipole interaction (DDI), which is long-range 
and anisotropic, plays a significant or even dominant role 
when compared to the short-range isotropic interaction. 
The special features of the DDI lead to fundamentally 
new physics in Bose-Einstein condensates [13, de- 
generated Fermi gases [HI, strongly-correlated atomic 
systems [3, quantum computation [l3], and ultra cold 
chemistry [l8|. Recently, time-of- flight experiments with 
Chromium condensates have allowed for the flrst obser- 
vation ever of dipolar effects in quantum gases [l9| . 

To a very good approximation the combined effects 
of the short-range interactions and the DDI can be un- 
derstood by means of a pseudopotential theory, which 
includes a contact interaction, characterized by a scat- 
tering length a, and the DDI itself [isl, |2o|. However, 
the correct value of a is in general not that in absence of 



DDI, but the result of the scattering problem including 
both the short-range and the DDI potentials. Hence the 
DDI may affect the value of a, even very severely in the 
vicinity of the so-called shape resonances [l^, [2^ [2lj . In 
addition, the combination of the DDI and the dressing 
of rotational excitations with static and microwave flelds 
may allow for the engineering of novel types of interaction 
potentials for polar molecules in 2D geometries [j^ . 

In this Letter, we analyze the scattering properties of 
dipolar gases in quasi-lD geometries. By solving the cor- 
responding scattering problem, including the short-range 
interaction, the DDI, and the trap potential, we obtain 
an effective ID pseudopotential consisting of a contact 
interaction with an effective ID coupling constant and 
a regularized dipolar potential. Similar as for the case 
without DDI we observe the appearance of a CIR, how- 
ever the position of the CIR is crucially modifled by the 
DDI. In particular, even for large values of a for which in 
a 3D environment the short-range interaction fully dom- 
inates the DDI, the latter can dramatically modify the 
ID scattering. As a consequence, the properties of the 
ID gas may be crucially modifled, leading to observable 
effects, which we discuss at the end of this Letter. 

In the following, we consider a dipolar gas of particles 
with (electric or magnetic) dipole moment c? in a quasi- 
lD geometry along the axial x-direction. The transver- 
sal 2/z-conflnement is given by an harmonic potential of 
frequency uj, whereas no trapping is assumed in the x- 
direction. We are interested in the scattering of two dipo- 
lar particles under these conditions, and hence (due to 
the separability of the Hamiltonian) we consider the cor- 
responding scattering problem in the relative coordinate 
r, given by the Schrodinger equation: 



(1) 



where = + , fi = m/2 is the relative mass, Vsh{f) 
is the short-range potential, and E = fku + h'^k'^ /2iJ., with 
k the axial momentum. Assuming the dipoles oriented 
along X, the DDI is given by: 
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and hence the DDI is attractive along x. This config- 
uration is particularly convenient for the analysis since 
it maintains the cylindrical symmetry of the problem. 
This symmetry may be also preserved if the dipoles while 
forming an angle (j) with the x-axis are rotated in the 
plane perpendicular to the trap axis fast enough to lead 
to an effective time averaged DDI [2^, similar to that 
of Eq. ([2]) but in which <P transforms into a<P , where a 
may range between 1 (0 = 0) and —1/2 (</> — 7r/2). Note 
that for a < the DDI becomes repulsive along the axis. 

If the system is assumed in the single-mode approxi- 
mation (SMA), i.e. only the transversal ground-state is 
considered, the wavefunction may be split as = 
f{x)(j)oip), where 0o(p) = expl-p"^ /2P]/^l, with P = 
h/{iJL<jj). After integrating over the transversal coordi- 
nates, the DDI leads to a delta-like interaction plus an 
effective ID DDI of the form: 
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where t — [x/Vf', crfc denotes the complementary error 
function. For distances x ^ I, V\d ex. 1/x^ as expected. 
Note that although the DDI is long-ranged and divergent 
at the origin, Vid{x) is regularized at a; = as well as in 
momentum space, being finite for fc = 0, leading us to 
use ([3]) in the pseudopotential below. 

To formulate an effective theory for the low-energy 
scattering of the dipolar particles a proper pseudopoten- 
tial must be introduced, which as shown below may be 
chosen as a combination of a regularized ID DDI and an 
appropriate contact potential. The motivation for this 
choice may be understood by considering, still within the 
SMA, the scattering of two dipoles at a short-range po- 
tential of range b. For d = 0, the psedupotential is given 
by giDd{x), where gm is a given coupling constant. How- 
ever, for d 0, the DDI dominates the asymptotic be- 
havior of the wavefunction for x >> b, which for fc = 
is given by the equation cPi/j/dx'^ ~ —{2apd'^/h'^)ip/\x\'^, 
whose general solution (for a ~ 1) is: 



^(x) = [AM2^2ld/x) + BYi{2^2ld/x)^ (4) 

where A and B are constants, Ji and Yi are Bessel 
functions, and Id = d'^p/h'^ is a length scale asso- 
ciated to the DDI. Asymptotically, ip{x 00) = 
A' [x/ Id + 2log{x/ Id)] + B', where A' and B' are con- 
stants. Note that the asymptotic form of the wavefunc- 
tion is not linear, and therefore the true potential cannot 
be replaced solely by a contact pseudopotential, as for 
d = 2J| . We choose thus a convenient pseudopotential 



providing the correct asymptotic behavior, formed by a 
contact interaction, and the regularized DDI ([3]): 



^e// = 9id6{x) + Vld{x) 



(5) 



where guj is obtained for a particular short-range poten- 
tial, after matching the logarithmic derivative at x = b. 



The SMA may become insufficient to describe the two- 
body scattering problem in a quasi-lD geometry. In 
particular, as discussed above, an analysis beyond the 
SMA [1, shows that for short-range interactions the 
scattering process in the presence of a transverse har- 
monic confinement potential undergoes a CIR, at which 
the coupling strength of the effective ID pseudopotential 
diverges. For a 3D contact pseudopotential ^^^-^^(^), 
the effective ID coupling constant is given by [§|: 
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and hence gio diverges when a ~ 0.68L It has been 
verified numerically that the shape and the position of 
the CIR does not depend on the details of the interatomic 
interaction, when the range of the interaction is much 
smaller than the radial confinement length 

In the following we analyze the effects of the dipo- 
lar interaction in the CIR by considering the scatter- 
ing process beyond the SMA. We solve the 3D scatter- 
ing (simplified by the cylindrical symmetry of the prob- 
lem) numerically. For simplicity, we consider a sim- 
plified model for the short range potential, namely a 
finite-depth potential well: Vsh = —VoO{ro — r), for 
which the 3D scattering length a is analytically known 
a = To [1 — taJl{^y mVoro / h) / y/mVoro / fi\ [25|. denotes 
the range of the interactions, which is kept small com- 
pared to the radial confinement length I (in the calcu- 
lations below we use tq — O.ll). In the following we 
assume that for the DDI considered the system is suf- 
ficiently far away from shape resonances, and hence we 
maintain the analytical value of a obtained in absence of 
the DDI. The 3D dipolar interaction is considered of the 



form 
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(r — ro), where we have imposed 

a cut-off at short distances to avoid divergences. This 
cut-off is physically justified since at short distances the 
short-range potential dominates. 

To evaluate the effective ID coupling constant guj 
we solve the 3D Schrodinger equation Q for a momen- 
tum fc = 0. Employing the cylindrical symmetry of the 
problem, we consider a 2D numerical grid for the ra- 
dial and axial coordinates, choosing the maximal radial 
value such that the wave function vanishes at the border, 
whereas the axial box has a length L » I » ld» r^. 
Note that the latter restricts our calculation to suffi- 
ciently small values of dipole moments to satisfy the 
ID condition I > Id- In order to increase the preci- 
sion near the scattering center, we employ a non-uniform 
grid 7y(i) =sech(A^i), where rj = p,x, and are prop- 
erly chosen. Since only the even parity wave function 
contributes to the low-energy scattering problem, we im- 
pose the boundary condition ^^g"'^-* = 0. For L ^ Id, 
the wave function can be written in a product form: 
ip{x,p) = f{x)<po{p). Imposing at x = L the logarith- 
mic derivative in a fully ID calculation, we evolve 
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FIG. 1: giD in units of h? / as a function of a/l for various 
values of d? /l^huj for (a) a = —1/2 (repulsive dipoles), and 
(b) Of = 1 (attractive dipoles). The bold line represents the 
analytical results of Ref. 
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FIG. 2: giB in units of / fil as a function of in units of 
I'^hcj for a fixed 3D scattering length a — 10.1/, for a = —1/2 
(dotted line) and q = 1 (solid line). 




FIG. 3: Quasi particle excitation energy E{k) in units of hu} 
as a function of kl, for 2nd? /l^ — Q.^huo, for ngio/huj = —0.4 
(solid), -0.43 (dashed), and -0.447 (dotted-dashed). 



the ID wavefunction from a; = L to a; = using ^1^(2;). 
From the logarithmic derivative of the ID wavefunction 
at X = we obtain giD- 

Fig. [T] shows the value of giu for different dipolar 
strengths (P/l^huj, as a function of the 3D scattering 
length. For vanishing dipole moment, our numerical 
results are in good agreement with the analytical for- 
mula ^ [8]. For sufhciently large values of a, the effec- 
tive coupling constant approaches an asymptotic univer- 
sal negative value — ^'^JJ^ ■ The DDI significantly modi- 
fies the behavior of gm as a function of a. With grow- 
ing d the position of the CIR shifts towards larger posi- 
tive values of a. At some particular value of the dipole 
strength {cP/l^ ^ O.Shuj), the position of the CIR is dis- 
placed towards a — +00, and as a consequence the CIR 
dissappears and gm monotonically increases with a. Fur- 
ther increasing \d\ shifts the CIR to negative values of a. 
When increasing d even further the CIR is found again 
for a > 0, scanning again all values of a until +00 and 
back from —00 to 0, and so on in a cyclic way. 

Note that in SMA a sufficiently large repulsive DDI 
on the X-axis would lead eventually to a shielding of the 
short-range potential, and gio would become an univer- 
sal function of d. However, in the calculations presented 
in this paper we consider the case in which I ^ Id- There- 
fore the anisotropy of the DDI becomes important, and 
the relative wavefunction can surround the repulsive axial 
barrier, and explore the short-range part of the potential. 

Note that in a 3D environment when |a| is sufficiently 



large the DDI becomes irrelevant compared to the short- 
range interactions. Remarkably, this apparently intu- 
itively obvious fact is not what occurs in a waveguide 
geometry. Fig. [2] shows the asymptotic value of guj for 
large |a|, which interestingly, due to the DDI can ap- 
proach a positive value, i.e. the gas acquires a repulsive 
character. This qualitatively differs from the behavior 
for d — 0, for which the effective ID contact interactions 
acquire for sufficiently large \a\ an universal attractive 
character [8|]. Hence, even if the short-range interactions 
fully dominate the physics in a 3D geometry, the DDI can 
dramatically change the properties of the dipolar gas in 
quasi-lD geometries, eventually changing the sign of gio- 
The prediction of this remarkable property introduced by 
the DDI (which should be experimentally observable for 
quasi-lD dipolar gases at a Feshbach resonance) can be 
considered the main result of this Letter. 

Once the pseudopotential is known, one writes down 
the many-body Hamiltonian for quasi-lD dipolar bosons 



H = dx-ii^x) 



where 4'(x) (^'^(x)) is the creation (anihilation) operator 
of bosons at x. At the CIR a Tonks gas with additional 
DDI would lead to a super- Tonks gas, with Luttinger 
parameter K < I [2^. For small values of |5i_d|, and al- 
though strictly condensation is prevented in ID, a finite- 
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size system at a sufficiently low temperature allows for a 
quasi-condensate [27]. This quasi-lD dipolar BEC may 
present a remarkable physics, as it becomes clear from an 
analysis of the dispersion law E{k) for axial excitations of 
momentum k on top of an homogeneous quasi-lD BEC: 



E{k) = ^e{k) [e(fc) + 2 (gii, - l/id(fc)) 



(8) 



where e{k) = h^k^/2m, Vid{k) = ^ [1 ^ ae'^T{0,a)] 
is the Fourier transform of Vidix), with a = k'^P/A. 
Note that |Vid(fc)| monotonically decreases with k, with 
l^id(0)| = 4:\a\(P/P. Hence, phonon stability demands 
giD > Vi£)(0), since otherwise the homogeneous quasi- ID 
BEC becomes unstable against the formation of bright 
solitons. However, even for guj > Vid{0), the homoge- 
neous BEC can become eventually unstable. This occurs 
when giD < 0, and q < (repulsive DDI along the axis). 
In that case the function guj — Vid{k) changes its sign for 
a sufficiently large fc, leading to the possibility of achiev- 
ing a roton (see Fi g. [3t in the spectrum at intermediate 
values of k 



H, 2^ 30|. When |5i_d| increases, the roton 
becomes deeper, and for a critical giD = ffc, the roton 
minimum reaches zero energy at a finite momentum kc, 
leading to an instability of the homogeneous quasi- ID 
BEC. For a = —1/2, and a given d, kc may be obtained 
from the equation (1 -I- ac)e'^^r[0, ac] = I + Phu/{4:(Pn), 
where CTc = k^P/A, whereas —2gcn/hw = (1 — ctc + 
a2e'^^r[0, (Te])/(1 - (1 + a,)e-^r[0, a,]). 

In summary, the DDI may play a crucial role in the 
physics of quasi- ID condensates, even for situations for 
which the short range interactions overwhelm the DDI in 
a 3D environment. The properties of the quasi-lD dipo- 
lar BEC crucially depend on the value of gio , which in 
turn will depend in a non-trivial way on the dipole and 
dipole-orientation due to the dipole- induced modification 
of the CIR. In particular, the DDI may change the sign 
of giD for a large 3D scattering length, changing com- 
pletely the physics of the quasi- ID dipolar condensates. 
For polar molecules the dipole strength may be modified 
by controlling the orienting electric field, and the DDI 
can be scanned from zero to large values. The combina- 
tion of this control and the modification of the induced 
CIR, should allow to scan close to the CIR all physics 
ranging from an homogeneous stable quasi-condensate, a 
bright-solitonic solution, a BEC with (eventually unsta- 
ble) roton minimum, and even a super- Tonks regime. 
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